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Abstract 

In this paper, we consider a drift-diffusion system describing the corrosion of an 
iron based alloy in a nuclear waste repository. In comparison with the classical drift- 
diffusion system arising in the modeling of semiconductor devices, the originality of 
the corrosion model lies in the boundary conditions which are of Robin type and 
induce an additional coupling between the equations. We prove the existence of a 
weak solution by passing to the limit on a sequence of approximate solutions given by 
a semi-discretization in time. 



1 Introduction 

1.1 General framework of the study 

In tiiis paper, we consider a system of partial differential equations describing the corrosion 
of an iron based alloy in a nuclear waste repository. At the request of the French nuclear 
waste management agency ANDRA, investigations are conducted to evaluate the long- 
term safety assessment of the geological repository of high-level radioactive waste. The 
concept for the storage under study is the following: the waste shall be confined in a glass 
matrix, placed into cylindrical steel canisters and stored in a claystone layer at a depth of 
several hundred of meters. The long-term safety assessment of the geological repository 
has to take into account the degradation of the carbon steel used for the waste overpacks, 
which is mainly caused by generalized corrosion processes. 

In this framework, the Diffusion Poisson Coupled Model (DPCM) has been developed 
by Bataillon et al [2] in order to describe the corrosion processes at the surface at the steel 
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canisters. It assumes that the metal is covered by a dense oxide layer which is in contact 
with the claystones. The model describes the development of the dense oxide layer. In 
most cases, the thickness of the oxide layer ranges from nanometers to micrometers and 
is much smaller than the size of the steel canister, which justifies a ID-modelling. The 
model is made of electromigration-diffusion equations for the transport of charge carriers 
(three species are considered: electrons, Fe^~^ cations and oxygen vacancies), coupled with 
a Poisson equation for the electrostatic potential. The interaction between the oxide layer 
and the adjacent layers (metal and claystone) are described by Butler- Volmer laws for 
the electrochemical surface reactions and potential drops. The model includes moving 
boundary equations. 

Numerical methods for the approximation of the DPCM model have been designed and 
studied by Bataillon et al in [Tj. Numerical experiments with real-life data shows the 
ability of model to reproduce expected physical behaviors. However, theoretical questions 
like existence of solutions and long-time behavior of the DPCM model have not been 
investigated yet. In this paper, we consider a simplified version of the DPCM model, 
already introduced in [T]. In this case, the interfaces of the oxide layer are assumed to be 
fixed and only two charge carriers are taken into account: electrons and cations Fe^"*". We 
will prove the existence of a global in time solution for this simplified model. 

1.2 Presentation of the model 

We recall here the simplified version of the DPCM model introduced in [T]. It is a di- 
mensionless model. The unknowns are the densities of electrons and cations, which are 
respectively denoted by N and P and the electrical potential denoted by ^ . The current 
densities are denoted by Jn and Jp; they contain both a drift part and a diffusion part. 
As we do not take into account the moving boundary equations in the simplified model, 
the domain describing the oxyde layer is the interval (0, 1). The equations describing the 
electrochemical behavior of the dense oxyde layer write, for t > 0: 

dtP + d^Jp = Q, Jp = -d^P-3Pd^^, in (0,1), (la) 

£dtN + d^JN = 0, Jn = -dxN + Nd^'^, in (0,1), (lb) 

-X^d,,^ = 3P-N + phi, in (0,1), (Ic) 

where phi is the net charge density of the ionic species in the host lattice (we will assume 
that it is a constant in all the sequel) and and e are positive dimensionless parameters 
arising in the scaling. Let us just mention that e is the ratio of the mobility coefficients 
of cations and electrons, then e <C 1. 

As the equations for the carrier densities (|la|) and (|lbp have the same form, we will use 
the synthetical form : 

Sudtu + dxJu = 0, Ju = -dxU - Zuudx'i', in (0,1). (2) 

For u = P,N, the charge numbers of the carriers are respectively = 3,-1 and we 
respectively have = 1)£- 

Let us now focus on the boundary conditions. Charge carriers are created and consumed 
at both interfaces x = and x = 1. The kinetics of the electrochemical reactions at the 
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interfaces are described by Butler- Volmer laws at interfaces. It leads to Robin boundary 
conditions on N and P. As in [1], we assume that the boundary conditions for P and N 
have exactly the same form. Therefore, for u = P,N, they write: 

-Ju = rl{u, ^) on X = 0, (3a) 
Ju = ri{u,^,V) on x = I, (3b) 

where y is a given applied potential (we just consider here the potentiostatic case) and 
rJJ and are linear and monotonically increasing functions with respect to their first 
argument. More precisely, due to the electrochemical reactions at the interfaces, and 
r^, for u = P,N, write: 

rl{s,x)=f3^^ix)s-jl{x), (4a) 

rlis, X, V) = PiiV - x)s - ^i{V - x), (4b) 

where the functions (/3^)i=o,i, (7u)i=o,i are expressed with the interface kinetic coefficients 
(m^, A;^)j=o,i; the positive transfer coefficients (o^, 6^)i=o,i, the maximum occupancy for 
octahedral cations in the lattice and the electron density in the state of metal N"'^. 
They write, for u = P,N: 

Pl{x) = mte"""*-^ + A;le^"<^ ^ = 0, 1, (5a) 

^o^x) = m°n"e-""''^, ^i{x) = A;>"e""<^ (5b) 

The boundary conditions for the Poisson equation take into account that the metal and 
the solution can be charged because they are respectively electronic and ionic conductors. 
Such an accumulation of charges induces a field given by the Gauss law. These accumu- 
lations of charges depend on the voltage drop at the interface by the usual Helmholtz law 
which links the charge to the voltage drop through a capacitance. The parameters A^fg^"^ 
and A^j'^'^ are the voltage drop corresponding to no accumulation of charges respectively 
in the metal and in the solution. Finally, the boundary conditions for the electric potential 
write: 

I- - ao(9a'^ = A'^l'"', onx = 0, (6a) 

^ + aid,,^ = V - A^I-f ^ on x = I, (6b) 

where oq and ai are positive dimensionless parameters arising in the scaling. 

The system is supplemented with initial conditions, given in L°°(0, 1): 

N{0,x) = No{xy, P{0,x)=Po{x), xeiO,l). (7) 

In all the sequel, we will assume that the interface kinetic coefficients are given constants 
and satisfy: 

mlklmlkl>0, foiu = P,N. (8) 
The positive transfer coefficients satisfy: 

alblalbie[0,l], ioiu = P,N, (9) 

and we assume that 

3P- _ + pi,i = 0. (10) 

Indeed, in the applications, we have p^i = — 5, = 2 and N'^ = 1, so that the relation 
()10p is satisfied. Then, it is expected that the solution to the corrosion model verifies 
< P{t, x) < P^ and < N{t, x) < N"^ almost everywhere. These estimates will be 
proved in Theorem 11.1) they are crucial for the proof of existence of solutions to the 
corrosion model. 
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1.3 Main result and outline of the paper 



Let us first note that the system of equations (|la|) - (fTc|) is a simplified case (linear case) 
of the so-called drift-diffusion system. This model is currently used in the framework of 
semiconductors device modelling (see for instance [ISJ ED ED] ) or plasma physics (see [7]). 
In this context the drift-diffusion model has been widely studied, from the analytical as 
from the numerical point of view. Let us refer toP[IUl[IIl[I21[I31[Il[ISllIZl[ISllliEIlElI 
for results on existence and uniqueness of solutions (using some fixed points theorem) and 
asymptotic limits of the model (as quasi-neutral limit, zero-relaxation time limit, ...). 
From a numerical point of view, different methods have been proposed and studied (see 
for instance [3 H El El ES E3] ) • 

In all these papers, the considered boundary conditions are Dirichlet, Neumann or 
mixed Dirichlet/Neumann boundary conditions (corresponding to the ohmic contacts and 
the insulated boundary segments of the semiconductor device). Then, the originality of the 
corrosion model described in this paper lies in the boundary conditions ([3D, ^ which are 
of Robin type, and induce an additional coupling between the equations. The authors have 
shown in [B] the existence of solution for such a model in the stationary case. Moreover 
in [1], a numerical scheme for the model ([I])-© is given and analyzed. 

Let us now state the main result of this paper: the existence and boundedness of a weak 
solution to system ([ID-dTD. 

Theorem 1.1. Assume Q, i^Uij and that Nq,Po G L°°(0, 1) satisfy 

0<Po{x)<P'", 0<No{x)<N'^, for almost every X e (0,1). (11) 
Let us also assume that A^'q^'^ and A^^^'^ verify 

^ 'l + log (a%k%ao) ) < Avl/g^^ < ^ (i + bg (a%k%ao) ) , (12a) 
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- ^ (l + log (6]vm]vai)) < AM/f ^ ^ ^ + {bWpC^i)) ■ (12b) 

Then the problem admits a weak solution (P, A^, ^) which is defined by P, N € 

L°°([0,r] X (0,1)) nL2(0,r;/7i(0,l)), L2(0,r;i/i(0,l)) for all T > and: 

For u = N,P, Vv3 G L^{0, T; H^{0, 1)), 

— Eu / udfifdxdt + e / uo{x)(p{0, x)dx — / / {—dxU — Zuudx^) dx^dxdt 
Jo Jo Jo Jo Jo 

'{/3iiV - ^{t, l))u{t, 1) - jUv - ^(t, 1))) ip{t, 1) 

(/30(*(t, 0))n(t, 0) - 7°(*(t, 0))) vit, 0)] dt = 0, (13) 



T 







+ 



and yip £ L2(0,r;i7^(0, 1)), 

rT rl rT )2 



I I dx^dxi^dxdt- I — (y- ^(t,l) - A'ff' )(^(t,l)dt 
Jo Jo Jo "1 

+ / — (^(t,0) - A^|;"")v7(t,0)dt = / / - N + phi)vdxdt. (14) 
Jo cto Jo Jo 



4 



Moreover the electron and cation densities N and P satisfy for all t G [0, T] and x G (0, 1) 

< P(t, x) < P"*, < iV(t, x) < N"". (15) 

In order to prove this Theorem, we consider a semi-discretization in time for (m)-©, as 
done by A. Jiingel and I. Violet in [TO] for the quantum drift-diffusion equations. We first 
show the existence of a solution to the semi-discrete problem. Then, we obtain some a 
priori estimates on the solution to the semi-discrete problem. These estimates will ensure 
the compactness of a sequence of approximate solutions. Finally by passing to the limit 
on the time step in the semi-discrete scheme, we obtain the existence of a weak solution 
for the problem ([I])-© satisfying (fTCj) . 

The paper is organized as follows. In Section 2, we present the semi-discrete problem 
associated to (H])-© and we prove the existence of a solution to this problem. In Section 
3, we show some a priori estimates satisfied by the semi-discrete solution. Finally, in 
Section 4, we establish the convergence of a sequence of approximate solutions and passing 
to the limit in the semi-discrete scheme, we prove Theorem 11.11 

2 The semi discrete problem 

In this Section, we propose a discretization in time of the system ([T])-([7D. Therefore, we 
introduce a uniform subdivision {tk)o<k<K of the time interval [0, T], for some T > 0. The 
time step is denoted by At and we have t^ = kAt. 

We set 

po = Po,A^° = iVo. (16) 
Starting from this initial condition, for sdl < k < K, we consider the following system: 

- A^a^.xi/fc = 3pfc _ AT^ + p^,, in (0,1), (17) 

|i (^fc+i _ ^fc) + d^Jt = 0, 4 = -d^u^+^ - Zuu^'+^d^^^ in (0,1), for u = P,N, (18) 
supplemented with the following boundary conditions: 

- ao^x^*^ = A^l'", on X = 0, (19a) 
+ aid,j,^^ = V - A^-f ^ on x = 1, (19b) 

- 4 = Pl{m^)u^+^ - 7°(^*^0, on ^ = 0, for n = P, N, (20a) 
= fil{V - ^^)n'=+i - 7^(y - ^'=), on x = 1, for u = P,N. (20b) 

Let us first note that, for each value of k, the system PT|) - ([2U|) is decoupled. Indeed, 
knowing P*"' and A^'^ , we define as the solution to (HID, ((121). As (fTTD is a classical elliptic 
equation, if the right-hand side belong to L'^{0, 1), ^''^ will be defined as the solution to 
a variational formulation (see below) of p7|) . p9p . Therefore, G H^{0,1) and as 
-^^0,1) C C°([0, 1]) the boundary conditions have a sense. Then, knowing , we 
define u^+^ = P'=+^ or A''+^ as the solution to 1^, ([20]). The equation ([18]) on is 
a linear convection-diffusion equation. Using the Slotboom change of variables [Q |2D] will 
lead to a classical elliptic equation for which we can also define a variational formulation 
and the associated weak solution. 
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Theorem 2.1. Under the assumptions of Theorem the semi-discrete scheme (|16p - 
()20p admits a unique weak solution (P'', N'', '^'')o<k<K with ^''^ G i7^(0, 1) for allO <k < 
K and P^, G H^{0, 1) for aUl<k< K. 

Proof. The proof will be done by induction. Let /c > 0, we assume that P^, G L°°(0, 1) 
(a fortiori satisfied if P'^, A^*^ G H^{0, 1) C C''([0, 1])) are known. This property is satisfied 
for k = 0. Then, the variational formulation associated to the problem on ^''^ (|17|) . (|19|) 
writes: 

X^a^i^^if) =b'l,{ip)y^ e H\0,1), (21) 

with 

Jo Oil CiQ 

bUv) = —{V- AM/D V9(l) + — A^I/r^(0) + f'iSP" - iV^ + PM)¥'dx. 
ai ao Jo 

But, the functional space H^{0, 1) can be equipped either with its usual norm: 



II^IIhi = / i^xu) dx + u dx, 
Jo JO 

or with the following equivalent norm (see [B]): 

= / {d^ufdx + u{Of + u{l)'^. 



Indeed, using the compact injection from H^{0, 1) to C([0, 1]) and an adaptation of Poincare 
inequality, we obtain the existence of k > such that : 

< |||n|||j;^i < k\\u\\h-i Vn G H^{0, 1). (22) 

K/ 

Using the norm ||| ■ Wh^, it is immediately clear that aq,{-,-) is a bilinear, continuous and 
coercive form in H^{0,1). Moreover, 6^ is a continuous linear form in H^{0,1). Then, 
the Lax-Milgram Theorem ensures that the problem (|2ip admits a unique solution ^''^ G 
i?^(0, 1). Thanks to the embedding from iJ^ (0,1) toC°([0, 1]), we also have ^f'' G C°([0, 1]). 

Let us now focus on the problem on = or A^'^+^ (fT5|) . ([201) • We first use the 

classical Slotboom change of variables introduced for instance in |2Uj : tt'^+^ = e~^"* Ct^^- 
The current density rewrites: J^' = —e~^'^^'^dxCu~^^ ^^^^ can consider the variational 
formulation: 

a'i{Cu,v)=bt{v),y^eHH0,l), (23) 



with 



a^{Ku,(p) = XT / e ''''^\ufdx+ [ e d^HudxH^dx 
Jo Jo 

+/3i{V - ^'=(l))e-^"*'(i)K„(l)(^(l) + /30(^'=(0))e-^"*'W^,(0)(^(0), 
buiv) = ^ /^Vx + 7,l(V^-^'(l)Ml)+7°(^'(0))^(0). 



Ai 







As is continuous in [0,1], it is bounded and e is also bounded and its lower 

bound is strictly positive. The hypothesis dHI) also ensures that the functions /3^, 7^, 7^ 
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are strictly positive. Therefore is a bilinear, continuous and coercive form in H^{0, 1) 
and 6^ is a continuous linear form in H^{0, 1). Then, the Lax-Milgram Theorem answers 
that the problem ()23|) admits a unique solution denoted Cu^^ ^ -f^^(0;l) ^-iid n*^'^^ = 
g-z„vi/*^fc+i g ifi(0, 1) is the unique weak solution to (USD, (IM . (l20bl) . This ends the 
proof of Theorem 12.11 □ 

Thanks to the semi-discrete scheme (|16p -(|20 p and Theorem l2.H we can define an approx- 
imate solution to the corrosion model ([I])-©. This approximate solution is a piecewise 
constant function in time. It is denoted by {PAt^ ^^At^^ At) and defined by 

PAt(.t,x) = P\x), NAt{t,x) = N'^ix), ^At(.t,x) = 

Vx G (0,1), Vt G itk-i,tk], yi<k<K. (24) 

Our objective is to prove the existence of a solution to the corrosion model by passing 
to the limit on a sequence of approximate solutions (Pa*! ^At^ ^At)At when At tends to 0. 
Therefore, we need to establish a priori estimates which will yield the compactness of the 
sequence of approximate solutions. 



3 A priori estimates 

In this Section we establish the a priori estimates satisfied by the approximate solution 
iPAt,NAt,'^At)- L°°-estimates on NAt and Pa*, an ^^(0, T, F^(0, l))-estimate on ^At 
and L2(0,r,i7^(0, l))-estimates on NAt and Pa*- 

3.1 L°°-estimates on PAt, NAt 

Proposition 3.1. Let us assume the hypothesis of Theorem M . 1\ and that the time step At 
verifies At < r with 

T = AVin f^-,-^V (25) 

for all k £ {0, K} we have 

< P'^(x) < P™, < N''{x) < N"", Vx G (0, 1), (26) 

and therefore 

< PAt(t,x) < P™, < NAt{t,x) < N"", Vx G (0,l),Vt G (0,r]. (27) 

Proof. The proof will be done by induction. First, we note that ()26|) is satisfied for 
A; = thanks to (|16|) and (|11|) . Let us now assume that the property is satisfied for some 
k G {0, ...,K — 1} , we will establish that both inequalities in ()26|) remain true at step 
k + 1. Therefore, we adapt the technique used by A. Jiingel and Y. J. Peng in [17j - 

Let u = P, N. We first prove the positivity of u'^"'"^. Multiplying (fTBj) by (u'^+^)~ = 
min(n'^^^,0) and integrating over (0,1), we obtain, after an integration by parts: 



At 
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with 

fu = J^(0)n^+^(0)- + ^(u'+\0)-)^d,^'{0), 
But, using the induction hypothesis and the fact that znzp < 0, we have: 



2A2 



1 / . X _\ 2 

fc+1 







> 



2A2 

and (pSj) imphes 



We remark that the condition on At (|25p and the induction hypothesis ensure that 

e.-^(n'"-n*^(x))>0 Vxe(0,l). (30) 

Let us now study the sign of and Therefore, we introduce, as in [T], the following 
functions : 

^uix) = 7°(x) - n-/30(x) + n-^(x - AM/^) Vx E M, 
= liix) - n-/3i(x) - n-^(x - AM^D Vx G M. 

As shown in [I], ^2 ^^'^ nonpositive functions under the hypothesis (|12p . Using the 

boundary conditions ([20]) . we have: 

fu = ^""'^2^!^"^' (^^(^'(0)) - u^p^^i^Ho)) - 7°(^'(o))) + ti^Ho))u'+Hor 

= ^ ' {euiv - ^'=(1)) - u^f^iiv - ^'{1)) - - ^'{1))) 

As the functions (Z^^, 7u)i=o,i are nonnegative on M and the functions (^^)j=o,i are non- 
positive on M for u = N, P, we deduce /° < and < 0. Then, ^ and ^ imply that 
_ Q Qj-^ -|^^ g^^^ obtain the positivity of u^"^^. 

Let us now prove that u^~^^ < n"*. Therefore, we multiply (|18|) by (u'^^-'^ — n™)"*" = 
max((ii^'"'"^ — n'"),0) and we integrate over (0, 1). It yields: 

+ y" [d^{u''+^ - u'")+y dx + zu J d^{u''+^ -u'^)+u''+^d^^''dx 

+ J^(l)(tx'^+i(l) - u"^)- - 4{0){u''+\0) - n'")+ = 0. (31) 
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But, as d^{u''+^ - = ^d^ (^((n'=+i - + 2n'"(«^+i - integrating 

by parts and using (fT7|) . we get: 

Jo 

- I {(u'^^Ho) - + 2n'"(n'=+i(0) - a^^^'(O). (32) 

Again, using the induction hypothesis and the fact that z^zp < 0, we have 



_z, 
2A2 



& lo (^''''^^ ~ ^ 2n"^(n'=+i - 7^"^)+) (2p(P*^ - P"^) + zn{N^ - iV™)) fix 



zlu"^ '■I 



Combined with ([?T|) and ([52]) . it yields: 







+ f^{d^{u^+^-u^Yy dx<gi + gl (33) 



with 

gl = J„(0)(n*^+i(0) - 1^"^)+ + ^ ((u'=+^(0) - n"^)+2 + 2u^{u^+\<d) - n'")+) 5,^^(0), 



i = -Ml){u''+Hl) - u"Y - y ((«'+'(!) - ^^")+' + 2n"^(n'=+i(l) - tx™)+) 9,.^'=(1). 



9u 

As we did for and , we use the boundary conditions (fTn|) and ([20]) and the definition 
of the functions (^^)j=o,i in order to rewrite ^^'^ 9u '■ 

9i = (^'^'(0)^-^")^^ {eui^'m - u-fiii^H^)) - 7°(*'=(o))) 

+ eO(M/^-(0))(n'-+'(0)-n™)+, 

+eu{v--^\mu^+\i)-u^)+. 

It is now clear that 5^ < and g\ <Q and therefore, thanks to the hypothesis on the time 
step (I25I), the inequahty ^ implies (n''+^ - n™)+ = and n'^'+i < for all x G (0, 1). 
It ends the proof of Proposition 13.11 □ 
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3.2 L2(0,T,//i)-estimate on 

Proposition 3.2. Under the assumptions of Proposition \ there exists a constant C 
depending only on the data of the problem (but not on At ) such that: 

III*'IIIhi((o,i)) < c", yi< k<K, (34) 

and 

II^At||L2(0,T,/fi(0,l)) - CVt. (35) 

Proof Applying (glD with = we get a*(^''',^'') = 6|(^''')/A2. Setting r/ = 
min(l, 1/ao, 1/ai), we have 

a^{^\^'')>r]\l^''f. (36) 
Using Young and Cauchy-Schwarz inequahties, we get 



^i^K^'^)! < Iw^'f + ^ I'isiP" - pn + 



_ N"')ydx 



+ ^(y_ A^f")2 + — (A^n^- (37) 
2r] 2ri 



Then, using the L°°-estimate ([25]) . we deduce ([5^ from pSj) and ([57]) . The obtention of 
([55]) is straightforward. 

□ 

Remark 3.1. Let us note that, due to the compact embedding from H^{0, 1) to C([0, 1]), 
the estimate (|34|) imply that ^''^ (for all k) is uniformly bounded on the interval [0, 1] by a 
constant depending only on the data of the problem. 



3.3 L2(0, T, //i)-estimates on Pa*, A^At 

Note that using (|27|) . we easily obtain L^(0, T; L^(0, 1)) estimates for P^t and A^a* which 
imply the weak convergence of a sequence of approximate solutions in L^(0,T;L^(0, 1)). 
But since we have to pass to the limit also in the space derivatives of Pa* a-^d A^a* ) it is not 
sufficient and we need estimates in L^{0,T; H^{0,1)). They are given in Proposition 13. 3[ 

Proposition 3.3. Under the assumptions of Proposition there exists a constant C , 
depending only on the data of the problem (but not on At ) such that 

l|-PAt||L2(0,T;//i(0,l)) < ||A'At||L2(0,T;Hi(0,l)) < C. (38) 

Proof. Let us multiply (jl8p by u^'^^ and integrate over (0, 1). Since 



we obtain 
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But, integrating by parts the term / dxJuU^~^^dx and using (fT7|) . we get 

Jo 



2At 



+ ^ j\u^+^)\zp{P^ - P-) + zn{N^ - N^))dx <hl + hi (39) 



with 



(',,fc+l('n^^2 

hi = (e°(*'(0)) - n-/3°(^'=(0)) - 7°(*'(0))) + 7°(^'(0))n'=+^(0), 

The sign properties of the functions /3^,7^ for i = 0, 1, ensure that 

hi < ^'^{^'m^'^'m and hi < ^iiV - ^\l))u'+\l). 

Then, we can multiply ([39]) by At and sum it over k E {0, . . . ,K — 1}. Using the L°°- 
estimate dSSD, we get 

l|9.«At||L^(0,T;L^(0,l)) < V f\u')\x)dx + ^Yl^t f\u'+')\u^ - u')dx 

I Jq Jo 

+ Y:^t[jl{^''{o))u^^\o) + jl{v-^Hm''^H^)) (40) 

Thanks to (jlip . ()26p . the continuity of the functions (7^)4=0,1 and Remark 13. 11 we obtain 

\\dxUAt\\L'2{0,T;L'2(0,l)) < C*) 

where C depends only on the data of the problem and not on At. As the L^{0, T; L^{0, 1))- 
estimate of u^t is a straightforward consequence of the L°°-estimate , it concludes the 
proof of Proposition 13. 3[ □ 



4 Proof of Theorem 11.11 

In this Section, we first establish some compactness results and the convergence of a 
sequence of approximate solutions (A^a*) -Pa* 1 ^Af)At when At tends to 0. Then, by passing 
to the limit in the scheme (|16p - (|2U|) . we prove that the triple (A^, P, ^') obtained as the 
limit of the sequence (A^t, -Pa*) ^At)At is a weak solution to ([I])-©. 



4.1 Compactness results 

The estimates proven in Propositions l3.2l and l3.3l implv the weak convergence of a sequence 
of approximate solutions (AAt)At) (-PAt)At and (^'At)At in -^^(0, T; //^(0, 1)), when At 
tends to zero. But, the weak convergence of (AAt)At and (^At)At in L^{0,T, {0, 1)) is 



11 



not sufficient to pass to the limit in the convective terms Jq Jq N^tdx'^ At^dx or in the 
boundary terms. We need some strong convergence. To prove this result, we will use 
a compactness argument due to Simon [23 ■ It is based on some estimates on the time 
translates of the approximate solutions. 

Let us introduce the shift operator. For ti = P, or we set 

(o-AfWAt) {x, t) = UAtix, t + At). 
It means that {(Tm'Um) {x,t) = u^~^^{x) for x G (0, 1) and t G (ifc-ii^fc]- 

Proposition 4.1. Under the assumptions of Proposition HOI there exists a constant C 
depending only on the data of the problem (and not on At) such that: 

\WAtPAt - -PAt|lL2(o,T;L2(o,l)) < CAt, \\cr At^At - ^At |Il2(o,T;L2(0,1)) < CAt. (41) 

and ||(TAt^At - ^AtllL2{o,T;L2{o,i)) < CAt. (42) 

Proof. Using equation ()18p . Minkowsky inequality and Holder inequality, we obtain for 
u = N, P 

\WAtUAt - UAt\\L^{Q,T;L^{0,l)) < — (\\9xCr AtUAt\\L2{0,T;L2{0,l)) 

+ |2;„|||9a:^'At||L2(o,T;L2(0,l))II^At'WAt||L°°(0,T;L°°(0,l)) 
At 



< — UI^At|lL2(o,T;//i(0,l)) 

+ N«|||^At|lL2(0,T;Hi{0,l))ll'"At|lL°°(0,T;L°°(0,l))) • (43) 

Using estimates ([27]), ([35]) and ([MD in dS]), we obtain dH]). 

By definition (dAt^At - ^At) {t,x) = ^''+\x) - ^^x) for t G {tk-i,tk]. Using (HZD, 
([T2D, we have 

- A^^^^ (^^=+1 - ^^•) = 3 (P'^^^ - P'^) - (iV'=+^ - A^'^) , xG(0,l), (44) 

(^fc+i _ ^fc) (0) - aodx (^^=+1 - ^'^•) (0) = 0, (45) 

Multiplying (|44|) by ^''^+^ — vl/'^ and integrating over (0, 1), we obtain, after an integration 
by parts: 



A2 



= sy - P'^) ("f - dx + J (iV^+^ - iV^') (^''^+1 - dx. (47) 

First, we note that the left hand side in (jlT]) is lower bounded by — with 

r] = min(l, 1/ao, 1/ai) and k defined in f21\i . Then, using Cauchy-Schwarz and Young 
inequalities, we obtain the following upper bound for the right hand side in ()47p : 
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It yields 

Summing over 1 < k < K and using the estimates (|41|) . we deduce (|42|) . It concludes the 
proof of Proposition I4.1[ 

□ 

Finally, thanks to the a priori estimates (Propositions IXTl 13. 2l and l3.3|) and the estimates 
on the time translates (Propositions 14.11 ) , we obtain the convergence of a subsequence of 
approximate solutions. 

Proposition 4.2. Under the assumptions of Theorem we consider a sequence of 
approximate solutions {PAti^At,'^At)At defined by the semi-discrete scheme (fTU|) - (pn|) 
and (IMI). Then, there exists {P,N) G L°°([0,T] x (0,1)) n L'^{0,T, H^{0,1)) and ^' G 
L^{0,T,H^ (0,1)) such that, up to a subsequence, (u = P, N) 

• "^At-^"^ strongly in L'^{0,T;C{[0, 1])), 

• UAt u strongly in L^(0, T;C([0, 1])), 

• dxUAt u weakly in L^(0, T; L^(0, 1)), (TAtUAt u weakly in L^(0, T; L^(0, 1)), 

• 5^^At ^ 9^'^ weakly in L'^{0,T; L'^{0, 1)). 

Proof. Firstly, due to Propositions 13.21 and 13. 3[ the sequences (^'At)At and {uAt)At are 
bounded in L'^{0,T; H^{0,1)). Then, there exist ^' and u in L'^{0,T- H^{0,1)) such 
that "if At tends to ^' and UAt tends to u in L^(0, T; L^(0, 1)). Therefore, up to a sub- 
sequence, we obtain the weak convergence of {dx'^At)At to dx"^ and of {dxUAtjAt to dxU 
in L2(0,r;L2(0,l)). 

Secondly, to obtain the strong convergence in L^(0, T; C([0, 1])) of sequences {^AtjAt 
and {uAt)At-, we use a compactness argument. More precisely, we use Theorem 5 given by 
J. Simon in [25]. Indeed, we have: 

i/i(0,l)cC([0, 1])cL2(0,1) 

with a compact injection from i/^(0, 1) to C([0, 1]) (due to the one dimension). Moreover 
due to Propositions 13.21 13.31 and 14.11 the sequence {fAt)At (for / = \I', u) is bounded in 
L'^{0,T;H^{0,1)) and 

\WAtfAt - /At|lL2(o,r;L2(o,i)) — > uniformly when At 0. 

Then the sequence {fAt)At is relatively compact in L^(0, T; C([0, 1])), which gives the 
strong convergence of {fAt)At in L^(0,r;C([0, 1])). 

Finally let us note that due to proposition 13.11 sequences {PAt)At and {NAt)At are 
bounded in L°°{[0, T] x (0, 1)) and we obtain their ^-weak convergence in L°^([0, T] x (0, 1)). 
Then their limits belong to L°°([0,T] x (0, 1)). 

This completes the proof of proposition 14.21 

□ 
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4.2 Passage to the limit 

In order to conclude the proof of Theorem II .![ it remains to prove that (P, N, ^) obtained 
as a hmit of the semi-discrete scheme (see Proposition 14. '2p satisfies and p4|) . 

Let us first recall the results of the previous sections. Theorem 12.11 gives the existence 
of a unique solution to the problem {u = N, P) 

{cTAtUAt - UAt) + dxJuAt = 0' JuAt = -dxicTAtUAt) - Zu{(JAtUAt)dx^At, (48) 

- X^d^^^At = 3PAt - N^t + Phi, (49) 

- JuAt = l3ui'^At)c7AtUAt - 7n(^At), X = 0, (50) 
JuAt = l^ii^ - '^At)(^AtUAt -ll{V- At), X = l, (51) 

^Ai - aod.^At = A^I-r', x = 0, (52) 
*Ai + aid^-^At = V- A^I-r, x = l. (53) 
The variational formulation of (IMI)-(IS3I) writes : for all ip G L'^{0,T; H^{0, 1)), 

f f CTAtUAt ^^t y^^^^j- _ f I (-d^aAtUAt - ZuGAtUAtdx^ At) d^^pdxdt 
Jo Jo Jo Jo 

- '^At{l))^AtUAt{l) - li{V - ^At{l))) V{1) 







+ 

+ (/3°(^At(0))aAtnAt(0) - 7°(^At(0))) ^{^)\ dt = 0, for u = N,P, (54) 



A2 



T 1 T \ 2 

/ / d.,^ Atdx^dxdt - [ —{V- ^At(l) - ^{'^)dt 
Jo JO Jo Ol 

+ / — (^Ai(0)-A^'n'^(0)(it= / [ i3PAt-NAt + Phi)^dxdt. (55) 
JO oo Jo Jo 



Note that 

Jo At ^ Jo Jo ^ ' At 

1 

uo{x)ip{0, x)dx 



Then, using proposition 14.21 and passing to the limit At tends to zero in (j54|) -(j55|) (up to 
a subsequence) , we get: 



7' 

Jo 



udt^pdxdt — / uo{x)(p{0, x)dx — / / {—dxU — Zuudx^) dxi^dxdt 
Jo Jo Jo 

fil{V-^{l))u{l)-jl{V-^{l))) ^{l) 
+ (f3'^{^{0))u{0) - 7°(^(0))) <^(0)1 dt = 0, for u = N,P, 



+ 
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T 1 T \ 2 

Jo Jo Jo Oil 

+ r — (^(0) - A^g^^) ^{Q)dt = r [\3P -N + phi)^dxdt, 
Jo ao Jo Jo 

for all if e L2(0, T; i7^(0, 1)). It is the expected result. 

Furthermore, passing to the limit in (|26|) gives for all t £ [0, T] and all x G (0, 1) 

< P{t, x) < P™, < N{t, x) < N"\ 

This ends the proof of Theorem 11.11 

5 Conclusion 

In this paper, we prove the existence of a weak solution to a simplified corrosion model. 
The boundedness of the electronic and cationic densities is also established. Let us note 
that the dimensionless parameter e in (jlbp is in practice very small because it is the ratio of 
the mobility coefficients of cations and electrons. It could be set to in the model. Then, 
the question of existence of a solution to the resulting model is still an open problem. 
Future work will focus on this question. The full understanding of the corrosion models 
reduced to a fixed domain will then permit to deal with the DPCM model on a moving 
domain (see [T]). 
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